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Metaheuristics have been acknowledged as an effective solution for many 
difficult issues related to optimization. The metaheuristics, especially 
swarm’s intelligence and evolutionary computing algorithms, have gained 
popularity within a short time over the past two decades. Various 
metaheuristics algorithms are being introduced on an annual basis and 


applications that are more new are gradually being discovered. This paper 
presents a survey for the years 2011-2021 on multiple metaheuristics 
algorithms, particularly swarm and evolutionary algorithms, to identify a 
nonlinear block-oriented model called the Hammerstein model, mainly 
because such model has garnered much interest amidst researchers to 
identify nonlinear systems. Besides introducing a complete survey on the 
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Metaheuristics , DS various population-based algorithms to identify the Hammerstein model, this 
Nonlinear system identification paper also investigated some empirically verified actual process plants 
Optimization results. As such, this article serves as a guideline on the fundamentals of 


identifying nonlinear block-oriented models for new practitioners, apart 
from presenting a comprehensive summary of cutting-edge trends within the 
context of this topic area. 
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1. INTRODUCTION 

System identification is vital in estimating a model that can mimic the attributes of a system [1]. 
Since past years, the identification and estimation of parameters for both nonlinear and linear systems have 
been vastly investigated by a range of researchers. Since most practical systems are nonlinear in nature, it has 
become more popular in recent times. The block-oriented models within the nonlinear system have been 
reckoned and mostly applied for automatic control and system detection. These block-oriented models are 
composed of a range of interconnections between static nonlinear and linear dynamic components. The three 
simple and widely-applied nonlinear block-oriented models to model practical systems refer to Hammerstein 
models and Wiener models, as well as the combination of both models [2]. 

In a traditional method of system identification, the reason for applying gradient-based algorithms 
for estimating parameters in a model is because it minimizes the mean square error (MSE) for both model 
and system [3]. For instance, an iterative algorithm based on a gradient descent algorithm was employed to 
identify the Hammerstein model [4]. In identifying the Wiener model and H-W model, Wang et al. [5] 
proposed gradient descent based and least square based iterative algorithm. The most generally used 
algorithm amongst all adaptive-type algorithms refers to least means square (LMS) [6]. In most cases of 
system identification, the error surface is multimodal and demands convergence into a sub-optimal 
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alternative by using the standard strategy of gradient descent [7]. This limitation can be overcome by using 
varied stochastic search approaches, such as the swarm intelligence based and evolutionary computing-based 
optimisation algorithms. 

Swarm intelligence algorithms refer to an efficient distributed intelligence framework for solving 
optimisation intricacies, which has been inspired by the mutual behaviour displayed by colonies of social 
insects and societies of other animals. Some optimisation algorithms that derived from the swarming 
behaviour metaphor are ant colony optimisation [8], particle swarm optimisation (PSO) [9], bacterial 
foraging optimisation (BFO) [10] and bee colony optimisation [11]. Evolutionary computation is the generic 
word over several optimisation algorithms influenced by the Darwinian principles of the power of nature to 
develop well-adapted living beings to their surroundings. The fields of genetic algorithm (GA) [12], 
evolutionary programming [13], and genetic programming [14] fall under the evolutionary computing 
algorithms, also known as evolutionary algorithms. The differential evolution algorithm (DEA) is also widely 
applied to address ongoing issues of global optimisation [15] and appeared to be a quite accurate optimisation 
technique for many distinct tasks. Many studies have also reported using such a type of optimisation 
algorithm to identify both linear and nonlinear systems. 

The pie chart in Figure 1 presents an overview of some popular block-oriented models used to 
identify nonlinear systems through the implementation of various metaheuristics algorithms. The pie chart is 
composed of three fractions that depict their usage popularity in identifying nonlinear systems; Hammerstein 
model leading with 53% (blue), followed by Wiener model with 32% (red), and the H-W model with 15% 
usage (green). It appears that most research favoured the Hammerstein model to identify nonlinear systems 
with a range of metaheuristics algorithms. 
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Figure 1. Pie chart of block-oriented models popularity using metaheuristics algorithms 


The bar chart in Figure 2 displays the relationships between different types of metaheuristics 
algorithms and their frequency of use for parameter estimation using block-oriented models. The 
metaheuristics algorithms are GA, PSO, DEA, glow-worm swarm optimisation (GSO), ant colony 
optimisation (ACO), BFO, colliding bodies optimisation (CBO), and brain storm optimisation (BSO). In the 
case of GA, Hammerstein and H-W models were used eight times each, while seven times for the wiener 
model. In PSO, Hammerstein, Wiener, and H-W models were used 17, 14, and three times, respectively. As 
for DE, Hammerstein, Wiener, and H-W models were used four, nine, and two times, respectively. Other 
metaheuristics algorithms were quite unfamiliar in identifying block-oriented models, as depicted in the bar 
chart. The bar chart shows that among all the metaheuristics; GA, PSO, and DE were mostly applied for 
parameter estimation. Thus, an extensive analysis is elaborated in this paper using three distinct optimisation 
methods (GA, PSO, and DE) to effectively identify the Hammerstein model. 

In this paper, we give an overview of the nonlinear Hammerstein system identification using several 
metaheuristics. Since the traditional methods have various limitations such as the creation of multimodal 
error surface, complexity in cost function minimization and inability to converse in global optima, the 
metaheuristics algorithms could be a good alternative for applying in nonlinear Hammerstein system 
identification. Metaheuristics can balance exploration and exploitation, it requires very few control 
parameters which will help to reduce the computational complexity and it has also the ability to converge in 
global optima. On the other hand, the Hammerstein model identification has been a hot research field because 
it can be successfully applied for the identification of real experimental plants and processes. Therefore, this 
paper could be provided a good research direction for the new researchers in the identification of nonlinear 
Hammerstein systems. 

The remaining part of the paper is as; section 2 discusses the general diagram of the Hammerstein 
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model and the review of the use of evolutionary and swarm optimisation algorithms to identify the 
Hammerstein model. Empirical verification of varied plants and processes is given in section 3, while 
section 4 depicts future research trends in this area. Lastly, section 5 concludes this study. 
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Figure 2. Bar graph: block-oriented models identification algorithms based on popularity 


2. RELATED WORK 
Scientific or information-based modelling, commonly known as system identification, is a 


significant aspect in the science and engineering domain. Comprehensive system identification development 
based on a linear model is vast and the related theories are mature. Researchers worldwide have been always 
interested in identifying both nonlinear and dynamic systems. For any control system design, system 
identification is the main component, and it is challenging to identify nonlinear systems in the present time. 
The nonlinear block-oriented systems, namely Hammerstein, Wiener, and H-W models, are interesting 
research areas amongst the academic and industrial environments. Among the three the Hammerstein model 
is the popular one which we can see from the pie chart in Figure 1. This section presents an overview of the 
block-oriented nonlinear Hammerstein model along with its block diagrams. 

The linear function in the Hammerstein model is preceded by nonlinear block, as shown in Figure 3, 
where u(t) is training input, x(t) denotes nonlinear system output, y’(t) represents Hammerstein model 
output, and v(t) signifies noise. The linear/nonlinear parameters are modified using metaheuristics algorithm 
for Hammerstein identification system, mainly to reduce cost function, for instance, MSE and least means 


square error (LMSE). From Figure 3, the nonlinear static block output is: 


x(t) = flu) (1) 
and y(t) is labelled as y(t) due to disruption by noise v(t). The input-output correlation is: 
y'(t) = H(s)x@) + v(t) (2) 
v(t) 
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Figure 3. Block diagram of block-oriented nonlinear Hammerstein model 


Figure 4 shows the conceptual structure of nonlinear system identification based on the 
Hammerstein model. Error minimization was determined to assess its objective function using (3): 
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J= RVO -IO dt (3) 


where J is defined as the objective function, y'(t) and f(t) are the real and estimated outputs, whereas t, is 
sampling time. The optimisation tool (algorithm) to update the Hammerstein model at the next iteration uses the 
variance between real and estimated outputs. This tuning step is iteratively performed until the last iteration. 
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Figure 4. Algorithm tuning process using the Hammerstein model 


This section also explains the Hammerstein model that is widely used to identify real nonlinear 
plants and processes [16], [17], including Thermoelectric cooler [18], [19], amplified piezoelectric tube 
actuators [20], forecast of wind speed [21], giant magnetostrictive actuators (GMAs) [22], boiler superheated 
steam pressure [23], boost converter of DC/DC [24], intelligent pneumatic actuator (IPA) [25], hydraulic 
mini excavator [26], ultrasonic motor [27], proton exchange membrane fuel cell (PEMFC) [28], turntable 
servo system [29], valve Stiction [30] and electrically stimulated muscle [31]. The Hammerstein model can 
accurately model actuators to solve control issues and strengthen nonlinearity system identification. This 
section reviews the Hammerstein identification scheme using soft and evolutionary computing methods. The 
two Hammerstein identification schemes are parametric and nonparametric. Nonlinearity is polynomial in 
parametric recognition system but modelled with probabilistic methods for the non-parametric system. The 
methods of Hammerstein identification are stochastic gradient method [32], least squares method (LSM) 
[33], subspace method [34], iterative method [35], as well as swarm and evolutionary algorithms. 

Akramizadeh et al. [36] initially proposed the GA for the identification of the Hammerstein system. 
They used GA for determining the correct nonlinear function structure and parameters, as well as the number 
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of zeros and poles of the linear transfer function. Lai et al. [37], the piezoceramic actuator is identified by the 
Hammerstein-based model and the parameters are identified using GA. They were used the Hammerstein- 
based model for the piezoceramic actuator for describing the frequency dependency. In 2015, Li et al. [38] 
proposed hybrid GA for the identification of the Hammerstein model and shows that hybrid GA has strong 
efficacy and robustness, has a good recognition effect, and is a viable solution for solving the nonlinear 
recognition problem. Qian et al. [29] identify the turntable servo system by using the Hammerstein model. 
They used GA to optimize the global control variable by minimizing the cost function, but their limitation is 
that the system properties differ in different working states. A nonlinear dynamic process called Blast furnace 
(BF) in ironmaking is identified and controlled in [39] where they used the Hammerstein model for the 
prediction of most essential quality (MIQ) indices. They employed the GA-based nonlinear model predictive 
control (GA-NMPC) method to demonstrate the applicability of the Hammerstein model for control design. 
Another plant called proton exchange membrane fuel cell (PEMFC) is modelled in [28] based on the 
fractional Hammerstein model and a fuzzy GA is adopted to acquire the best fractional order of the 
Hammerstein system. Next, a nonlinear multi-input single-output (MISO) system identification method using 
Hammerstein model was presented in [40], which involved fractional transfer function and modified radial 
base function NN (MRBENN) as the linear dynamic component and a static nonlinear subsystem, in which 
modified GA defined the fractional orders, as well as the width and centre of MRBFNN, apart from 
identifying the initial prediction of other unidentified parameters. Genetic algorithms (GAs) are used to 
identify parameters in a nonlinear Hammerstein controlled autoregressive (NHCAR) system but the 
optimization process necessitates additional execution time, generation counts, and function counts [41]. A 
similar approach is used for identifying the continuous-time Hammerstein system in [42]. In the paper they 
hybrid the GA with the recursive least-squares (RLS) method and their proposed method shows a superior 
identification performance. In another study, the Hammerstein controlled auto regressive auto regressive 
moving average (HCARARMA) system is proposed and the parameter is estimated by using DE, GAs, 
pattern search (PS) and simulated annealing (SA) algorithms [43]. They show that the complexity measures 
for GAs are higher than the rest of the optimization mechanisms. Similarly, a nonlinear Hammerstein 
controlled auto regressive auto regressive (NHCARAR) system is proposed in [44] and the parameters are 
identified by using a backtracking search algorithm (BSA), DE and GAs. In the study, it was found that the 
BSA and DE algorithms outperformed the GAs, with BSA providing equivalent accuracy at the expense of 
greater complexity than DE. 

PSO and its variants are the most popular identification method for the Hammerstein system. Using 
PSO, parameter estimation, and model selection was performed for the Hammerstein model. A method to 
quantify valve stiction in control loops using PSO was reported in [30]. The PSO from input-output data were 
used to estimate the parameters of the Hammerstein model by decreasing the errors between real and 
identified model outputs. The immunized PSO, which exerted more CPU load, was found to enhance 
modelling precision more than PSO did, while the clonal PSO gave PSO-like modelling precision and CPU 
load. Enhanced PSO gave better output than conventional PSO did [45]. For iterations of a MISO 
Hammerstein system identification task, the improved PSO (IPSO) gave more rapid convergence and better 
accuracy than the conventional PSO did, while neglecting runtime iteration. A work comparable to [45] was 
presented in [46], with the distinction that in place of mutation, the IPSO contained dynamically reduced 
inertia. To depict nonlinear block, the Hammerstein model was applied with non-uniform rational NN B- 
spline (NURB) [47]. NN shaping parameters were estimated by PSO and the other parameters using the over- 
parameterization method [48], but the evolutionary algorithms used for training were unreported. Upon using 
both maximum likelihood adaptive PSO (ML-APSO) and maximum likelihood PSO (ML-PSO) for the 
Hammerstein model, ML-APSO gave better convergence speed and modelling precision than ML-PSO and 
recursive least square (RLS) algorithm did [49]. Although convergence speed and efficiency were achieved 
by RLS-PSO, no comparison was made with PSO and RLS [50]. The CBO algorithm, inspired by the 
concept of a collision between bodies [51], was implemented in [52] for Hammerstein system modelling. 
Similarly, the performance displayed by CBO in [53] was better in terms of CPU time and MSE when 
compared with outcomes retrieved from CSA [54]. The Hammerstein model, along with other systems, were 
identified using a modified fish swarm algorithm (MFSA) [55], whereby MFSA displayed superior efficiency 
to GA, FSA, and PSO algorithms, while dismissing the convergence features of the algorithms. A similar 
algorithm was used to model a nonlinear dynamic system [56]. The CSA-based Hammerstein model was 
applied with FLANN and adaptive IIR filter modelling nonlinear and linear blocks, respectively [54]. As a 
result, the performance exhibited by CSA was superior to DE and PSO considering CPU time and lower 
MSE. Next, the Hammerstein model with GA managed to further decrease CPU time and MSE, when 
compared to GA alone, DE, and PSO [57]. Upon comparing two models of Hammerstein trained using DE, 
CSA, RLS, and GA; MPN with finite impulse response (FIR) filter (MPN-FIR) and volterra second-order 
with FIR (SOV-FIR), DE offered the best outcomes [58]. Using the Hammerstein model, the PSO and 
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Weierstrass estimation theorems were applied to estimate the parameters and optimize the nonlinear system 
[59]. By comparing the modeling performance between the stochastic gradient algorithm and the modified 
method, PSO had improved the modeling accuracy over other comparative methods. A discrete-time 
fractional-order Hammerstein model is identified using the PSO algorithm [60]. In the presence of noise, 
their simulations confirm the efficiency of the method, but the result's accuracy decreases with a large 
amount of noise. Another proposed version of PSO is called the craziness based particle swarm optimization 
(CRPSO) algorithm for solving the local optima problem of PSO [61]. They apply their CRPSO algorithm 
for the identification of the discrete-time Hammerstein model. An adaptive particle swarm optimization 
algorithm (APSO) was proposed to identify the single input single output (SISO) Hammerstein system [62]. 
To improve the convergence speed of the conventional PSO, they applied an evolutionary state estimation 
technique. A new algorithm based on the mixed-integer particle swarm optimization (IR-PSO) for the 
estimation of the integer orders and the real coefficients of the block-oriented nonlinear Hammerstein 
autoregressive with exogenous input (NLHARX) system was proposed in [63]. Although PSO has gained 
attention from these researchers due to its searching capability, as well as its novelty, PSO can easily suffer 
from local optima stagnation when handling optimization problems with multiple local optima. Additionally, 
proper parameter settings will greatly influence the searching performance of PSO. Hence, the local optima 
entrapment and premature convergence of PSO often occur in numerous engineering optimization problems [64]. 

Xiong et al. initially introduced DE for the Hammerstein model to identify the position control in 
the ultrasonic motor [65]. The DE algorithm is used in this context to identify both the orders and parameters 
of the model. A study proposed the adaptive mutation DE (MDE) to estimate the parameters in the 
Hammerstein model. In contrast, another study optimized the Hammerstein model using DE to identify the 
system [66]. The MPN and the FIR models were used for nonlinear and linear elements, respectively. But 
they claim that, in terms of run times by DE, another algorithm has not always been very successful. This is a 
disadvantage of DE. Many studies considered that the modelling of polynomial nonlinearity is like that of a 
real one, which is uncertain in practical scenarios. Thus, determining the optimum extent of the polynomial 
nonlinearity is more practical. Next, the Hammerstein model with GA managed to further decrease the CPU 
time and MSE, when compared to GA alone, DE, and PSO [57]. Upon comparing the two models of the 
trained Hammerstein using DE, CSA, RLS, and GA, MPN with finite impulse response (FIR) filter (MPN- 
FIR) and volterra second-order with FIR (SOV-FIR), DE offered the best outcomes [58]. In the study, it was 
found that DE required fewer control parameters compared to a well-known method, particle swarm 
optimization (PSO). These control parameters included population size, step size, and crossover rate. 
Furthermore, the study made by Vesterstrom and Thomsen [67] claimed that DE is better than PSO and EAs 
as it can obtain the lowest fitness solution for most of the problems reported in the literature. It was again 
reported that the convergence speed of DE is significantly faster than GA [58]. Additionally, DE is robust as 
it can consistently reproduce the same solutions over many simulations compared to PSO. This nonlinear 
USM of the Hammerstein model was introduced to properly categorize a nonlinear motor [27]. To improve 
both efficiency and performance of the model, a combination of the DEA and Hammerstein model was 
prescribed. A nonlinear Hammerstein controlled autoregressive (NHCAR) system was proposed, where the 
DE and GA were used for parameter estimation [44]. De also used to identify the Hammerstein system in 
[68], where they represent a real plant called electrically stimulated muscle (ESM) by using the Hammerstein 
structure. They claim that the DE method outperforms other metaheuristics such as GAS, PSO, pattern search 
(PS) and simulated annealing (SA). Recently, another real plant (ultrasonic motor) is modelled with the 
Hammerstein system and the DE algorithm is used to identify the model parameters [69]. A timeline of 
studies (2011-2021) involving the Hammerstein model incorporating evolutionary computing and swarm 
algorithms is given in Figure 5. 

A few other metaheuristics algorithms also reported the identification of the Hammerstein system. 
The cuckoo search algorithm (CSA) based Hammerstein model was applied with FLANN and adaptive IIR 
filter modelling in nonlinear and linear blocks, respectively [54]. As a result, the performance exhibited by 
CSA was superior to DE and PSO considering the CPU time and lower MSE. The colliding bodies 
optimization (CBO) algorithm, inspired by the concept of a collision between bodies, was implemented for 
the Hammerstein system modelling [52]. The CBO gave better modelling precision and less CPU time than 
bacterial foraging optimization and adaptive PSO. Similarly, the performance displayed by CBO [53] was 
better in terms of CPU time and MSE when compared with outcomes retrieved from CSA [54], PSO, GA, 
and RLS, regardless of the colourful noise. A new meta-heuristic orthogonal colliding bodies optimization 
(OCBO) was proposed for the identification of the Hammerstein plant [70]. However, the drawback was the 
computational time of their proposed OCBO model is higher. Another metaheuristic called gravitational 
search algorithm (GSA) is proposed by Erik et al. in [71] for the identification of the ANFIS-Hammerstein 
model. They claim that the GSA algorithm is outperformed by the PSO and DE by avoiding critical flaws 
such as premature convergence to sub-optimal solutions. In another study, novel hybrid metaheuristics called 
average multi-verse optimizer and sine cosine algorithm (AMVO-SCA) is proposed for the identification of 
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the continuous-time Hammerstein model [72]. They showed that their proposed method provides good 
identification accuracy than several well-known metaheuristics such as PSO, GWO, MVO and SCA. 


> Hybridized GA and PSO 
> NURB with PSO 
> Modified Adaptive PSO (MAPSO) 
> Enhanced PSO 
> Maximum Likelihood PSO (ML-PSO) 
> Maximum Likelihood Adaptive PSO (ML- 
PSO), 
> Adaptive Mutation Differential Evolution 
Algorithm (MDE) 
> Recursive Least Squares PSO (RLS-PSO) 
> Colliding Bodies Optimization (CBO) 
> Weierstrass approximation with PSO 
> Gravitational Search algorithm (GSA) 
> Cuckoo Search Algorithm (CSA) 
> Craziness based PSO (CRPSO) 
> Differential Evolution Algorithm (DEA) 
> Colliding Bodies Optimization (CBO) 
> Orthogonal Colliding Bodies Optimization 
(OCBO) 
> MRBENN with GA 2017 
> NLHARX with PSO . 
> Gravitational Search Algorithm (GSA) 2018 
> Fuzzy Genetic Algorithm 
> Differential Evaluation Algorithm 
> Artificial Bee Colony algorithm (ABC) 
> Position control using DEA 
> Functional Genetic Algorithm 


Figure 5. Timeline: Hammerstein model identification algorithms 


3. EXPERIMENTAL WORKS ON IDENTIFYING HAMMERSTEIN MODEL 

Empirical studies on nonlinear modelling and the identification of a real nonlinear plant are briefly 
discussed in this section. The experimental verification of such plants and processes is displayed in Table 1. 
These process plants can be classified into several groups, such as chemical plants, mechanical plants, 
industrial plants and engineering plants such as amplified piezoelectric tube actuators [20], intelligent 
pneumatic actuator (IPA) [25], electrically stimulated muscle [31], ultrasonic motor [27], proton exchange 
membrane fuel cell (PEMFC) [28], turntable servo system [29], hydraulic mini excavator [26], twin rotor 
system (TRS) and flexible manipulator system (FMS) [72], and valve stiction [30]. 

Both the modelling and compensation of hysteresis nonlinearities of piezoelectric tube actuator 
using inverse rate-independent Prandtl-Ishlinskii model (RIPI) and inverse rate-dependent prandtl—ishlinskii 
model (RDPI), along with a controller [20]. For their identification, the Hammerstein model was applied. The 
IPA system that applied the Hammerstein model had been based on the RLS algorithm [25] and they use a 
discrete-time ARX model to represent the linear component of the Hammerstein model. The nonlinear 
Hammerstein model of the ultrasonic motor was reported in [27]. The PEMFC modelling that used the 
Hammerstein model [28] showed that the fractional Hammerstein PEMFC model with subspace 
identification algorithm not only avoided the analysis of internal complex mechanism but also accurately 
described the stack and nonlinear fractional characteristics. The identification of nonlinear predictive 
functional controller (NPFC) and nonlinear model, based on Hammerstein model for turntable servo system, 
was reported in [29]. A Hammerstein structure is used to represent the isometric response in the paper [31], 
which considers modelling of an electrically activated muscle, where they proposed the ARLS algorithm for 
identifying the parameters. Jui and Ahmad [72], a continuous-time Hammerstein model is developed based 


Bulletin of Electr Eng & Inf, Vol. 11, No. 1, February 2022: 454-465 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 Oo 461 


on AMVO-SCA and validated by modeling two real world plants such as TRS and FMS. In the Table 1, the 
objective function equation of the process plants with its parameter identification algorithm, number of 
design parameter from both nonlinear and linear subsystem of the Hammerstein model and the type of 
Hammerstein model also reported. From the tabular value it can clearly see that, most of the existing 
processes and plants identification algorithm identify an average number of design parameters and the 
majority of the previous study represent the model in discrete-time domain. 


Table 1. Experimental verification table of Hammerstein model identification 


Parameter Number Types of 
Plants and processes estimation of design yP f Objective function equation Year Ref. 
3 Hammerstein 
algorithm parameter 
Valve stiction PSO 8 Discrete time M 2011 [30 
Yodo - 900)" 
k=1 
Electrically stimulated Alternately 5 Discrete time 1 N 2012 [31] 
muscle RLS (ARLS) TT Or -7) 
t=1 
DC/DC boost converter N/A 8 Discrete time 5” ae e(o) 2015 [24] 
k=kj 
Amplified piezoelectric RLS 3 N/A N/A 2016 [20] 
actuators (APA) 
Turntable servo system CLPSO & 8 Discrete time ule 2016 [29] 
PSO > b@ -900P 
k=1 
Intelligent pneumatic RLS 8 Discrete time ly -9| 2017 [25] 
actuator (IPA) 1- y-y * 100 
Hydraulic mini excavator N/A N/A Discrete time sA apasih Jo 2017 [26] 
min zu Hü+guü+v 


Ultrasonic motor (USM) PSO and DE 12 N/A i 2019 [27] 


: YD = yi)]?/h /4 
k=1 


i=1 


Proton exchange FGA 4 Discrete time N. 2019 [28] 
membrane fuel cell min I; -= V;)? + (P; — P)?] 

(PEMFC) i=1 

Twin rotor system (TRS) AMVO-SCA 5 & 13 Continuous N: 2021 [72] 
and flexible manipulator time 3 (y(yts) — ¥(nts))? 

system (FMS) =0 


4. FUTURE RESEARCH TRENDS IN HAMMERSTEIN MODELS IDENTIFICATION 
This work surveyed several important metaheuristics algorithms for the identification of the 

Hammerstein model as they are described in different studies. Despite a limited theoretical foundation, the 

benefits of metaheuristics are vast. Nevertheless, several issues need to be highlighted to further determine 

their full potential. Metaheuristics assessments were mostly based on empirical comparisons. The following 
lists some significant related research issues that require further works in future. 

— Despite the abundance of the existing literature on the Hammerstein model for system identification, the 
vast majority has focused on the SISO system. Nonetheless, works related to a multiple-input multiple- 
output (MIMO) system is in scarcity. The real objects are always nonlinear in MIMO. In control 
practices, with the increasing requirements to deal with the complexities of modern technology, MIMO 
systems have more practical and wider applications than SISO systems do. 

— Fractional order models have an important interest owing to their ability to represent many physical 
phenomena and being useful for modelling electrodynamics processes. These models have a remarkable 
advantage by representing the system with a low number of parameters. Identification of fractional 
order systems seems to be an active area within the research domain and mostly places focus on linear 
cases. 

— Despite the presence of several algorithms in addressing several drawbacks, the future views 
hybridization as a powerful approach to combine the benefits of multiple algorithms. This hybridisation 
approach, nonetheless, is still in its trial-and-error phase although attempts can be noted as back as 
many years ago. 

— Itis worthwhile to investigate a few current metaheuristics algorithms that use a small number of 
coefficients to reduce the time spent determining the best setting for the coefficients. 


Metaheuristics algorithms to identify nonlinear Hammerstein model: a decade survey (Julakha Jahan Jui) 


462 o ISSN: 2302-9285 


— Most earlier studies reflect discrete-time models, whereas most real-world systems are expressed in the 
continuous-time domain. It is worth researching a continuous-time Hammerstein model utilizing the 
metaheuristics method [73], [74]. 


5. CONCLUSION 

The domain of system identification is related to control engineering that uses mathematical models 
of nonlinear systems based on input-output observations. Swarm intelligence and evolutionary computing are 
essential in this study area, whereby algorithms based on swarm intelligence have a role in nearly all science 
and engineering research areas. An extensive review is presented in this paper regarding the implementation 
of various metaheuristics algorithms for the identification of the Hammerstein model. This paper introduces 
the flow diagrams of common metaheuristics algorithms considering nonlinear Hammerstein system 
identification that may serve as guidance for those also in the artificial intelligence topic area. Some 
empirical studies on identifying the Hammerstein model are tabulated, including several uprising issues. It is 
hoped that more studies would investigate computations that involve nature and swarm intelligence to 
address a range of optimisation issues in actual implementations. 
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